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_. Abstract 

CN I We introduce a simple geometric model which describes the kinetics of fragmentation 

of (i-dimensional objects. In one dimension our model coincides with the random 
^ . scission model and show a simple scaling behavior in the long-time limit. For d > 

^ ■ 1, the volume of the fragments is characterized by a single scale 1/t, while other 

t^^ ' geometric properties such as the length are characterized by an infinite number of 

O ' 

■r:;j- ' length scales and thus exhibit multiscaling. 
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I. Introduction 

The phenomenon of fragmentation which occurs in numerous physical, chemical, 
and geological processes, has attracted a considerable recent interest. Fragmentation 
can be exemplified by polymer degradation, grinding of minerals, atomic collisions 
cascades, energy cascades in turbulence, multivalley structure of the phase space of 
disorder systems, etc. [|I], 0, ^ ^, ||, |^, |^. In general, fragmentation is a kinetic 
process with scattering, breaking, or splitting of particular material into smaller 
fragments. With such wide-ranging applications it is natural to try to abstract the 
essential features of fragmentation and to model them as simple as possible. One 
characteristic feature of these cascade processes is that fragments continue splitting 
independently. This allows one to describe the evolution by linear rate equations. 
Another restriction which is used in almost all studies of fragmentation is the implicit 
assumption that fragments may be described by a single variable, say their mass or 
size. The simplest model satisfying these restrictions is the so-called random scission 
model 1^, |, p. In this model, the distribution of sizes is described by the integro- 
differential equation. 



dP{x,t) 
di 



-xP{x,t) + 2jdyP{y,t), (1) 



where P{x,t) is the concentration of fragments of mass (size) x, x-mers, at time t. 
The loss term on the right-hand side represents the decrease of x-mers due to binary 
breakups. The probability of breaking at every point is assumed to be constant in the 
random scission model and hence the overall rate at which an x-mer breaks is equal 
to X. The gain term in Eq. (|1|) represents the increase of x-mers due to breakups of 
longer fragments. The general solution |^, ||, |[ to Eq. (|l]) is 

P{x,t) = e-^* Ipo{x) + J dyPo{y)[2t + t\y - x)]\ . (2) 

In the long-time limit, this exact solution approaches the scaling form 

oo 

P(x,t) = CtV^*, C = JdyyPo{y), (3) 



if we keep xt finite while taking a limit t —>■ oo and x — > 0. 

The random scission model is a representative example of "one-dimensional" frag- 
mentation processes in which the fragments are described by a single variable. The 
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kinetics of such fragmentation processes is now well understood and numerous ex- 
plicit and scaling solutions have been found @, ||, P, |^, |10|, |Tl], |12|, |13], Q . 



The geometry of fragments clearly influences the fragmentation processes. How- 
ever, it was ignored in so far studied models. In this paper, we introduce simple 
kinetic models describing the splitting of two dimensional and more generally d- 
dimensional objects. We find that multiscaling occurs for dimensions larger than 
one. The rest of this paper is organized as follows. In section II, we present the 
two-dimensional model and analyze the behavior of the moments of the size distri- 
bution of the fragments. Furthermore, we investigate the area distribution of the 
fragments an show that it exhibits ordinary scaling. In Section III, we generalize 
the asymptotic results to arbitrary dimensions and show that multiscaling occurs in 
higher dimensions as well. In Section IV, we introduce a two-dimensional isotropic 
fragmentation process. Numerical study of this process suggests that it belongs to a 
different universality class. 

II. Fragmentation in two dimensions 

In close analogy with the one-dimensional fragmentation process we study the 
following process in two dimensions. The fragmentation event takes place at arbitrary 
internal point of the rectangular and gives birth to four smaller rectangulars as 
illustrated in figure 1. The distribution function P{xi,X2',t) describing rectangulars 
of size Xi X X2, is governed by the following kinetic equation, 

^''''''''^'- = -x,x,P{x,t)+Af fdy,dy2P{yuy2;t). (4) 



dt 

XI X2 



Note that Eq. (^ implies the conservation of the total area, 

dxidx2XiX2P{xi, X2;t) = const. (5) 





To analyze Eq. (^) we introduce the double Mellin transform of the distribution 
function P{xi,X2;t), 

oo oo 

M{si,S2;t) = / dxidx2xl^~^X2^~^P{xi,X2;t). (6) 





Figure 1: The Fragmentation Process. 



The functions M{si, S2', t) at fixed si and S2 will be called the moments. By combining 
Eqs. @ and (^) we arrive at the equation 



dM{si,S2;t) 
dt 



4 



S1S2 



1 M(si + 1,52 + 1; t). 



(7) 



A surprising feature of Eq. (j^) is that it implies the existence of an infinite number 
of conservation laws. The moments M(si, S2', t) with si and S2 satisfying the relation 
S1S2 = 4 are independent of time. Thus in addition to the conservation of the total 
area there is an infinite amount of hidden conserved integrals. These integrals are 
in fact responsible for the absence of scaling solutions to Eq. (^). Indeed, with the 
scaling solution P{xi,X2',t) = t'^Q{t^xi,t^X2), implies an infinite amount of scaling 
relations, w = z{si + S2) at S1S2 = 4, which cannot be satisfied by two scaling 
exponents, w and z. 

We will solve Eq. (^ by Charlesby's method p[ (for more recent applications 
of this method see, e. g., |l^, |15]). For the random scission of the unit square, 
P{xi,X2;0) = (5(xi — l)5(x2 — 1), or equivalently M(si, S2; 0) = 1. By iterating Eq. (|^) 
one can compute all derivatives of M(si, S2', t) at t = and then find M(si, S2', t) from 
the Taylor's series, M(si, S2] t) = M(0) + tM'(O) + t'^M"{0)/2\ + t^M"'{0)/3\ + . . .. 
This gives a solution in terms of a generalized hypergeometric function ||16||, 



M(si, S2; t) = 2F2{a+, a_; si, S2; -t), 



(8) 



with 



a± = a±[si,S2) 



Sl + S2 



± 



Sl - S2 



+ 4. 



2 V V 2 

Computation of first few moments gives M(l, 1; t) = N{t) = 1 



(9) 



3t for the total 



number of fragments, N(t); M(2,2;t) = 1 for the total area; and M(3,3;t) = j + 
3^ + e~* (I ~ ;!; + 342) for the next diagonal moment. The first moment can be easily 
understood. The rate of creation of fragments is equal to 3 since every fragmentation 



event introduces 3 additional rectangulars, and hence, the total number of fragments 
is 1 + St. These results suggest the following power-law asymptotic behavior of the 
moments M(si, S2]t), 

M(si,S2;t)~A(si,S2)t-"(^^'^^^ (10) 

Substituting this asymptotic form into Eq. (0) yields the difference equations for the 
exponent a{si, S2), and for the prefactor A{si, S2), 

a(si,S2) + l = a(si + 1,S2 + 1), 

4 X (11) 



a{si, S2)A{si, S2) = 1 A{si + 1, S2 + 1). 

\ S1S2/ 

With the boundary conditions, a(si, S2) = and A{si, S2) = 1 at S1S2 = 4, Eq. (11) 
are readily solved to give 



a{si,S2) =a_(si,S2) = — v( \ J + ^' 

A(s,,s2) r(.i)r(.2)r(a+-a-) ^^^^ 

r(a+ - si)r(a+ - S2)r(a+) ■ 

The preceding formulas, Eqs. (10) and (12), might be established rigorously from the 
asymptotic behavior of the generalized hypergeometric functions. 

For ordinary scaling distributions the exponent a{si,S2), describing the asymp- 
totic decay of the moments is linear in the variable Si + S2- However, for two- 
dimensional fragmentation this exponent depends also on the variable si — S2. This 
manifests the non-trivial scaling properties of this process. One can also compare 
the average value of Xl^X2^ defined by 

00 00 

J J dxidx2Xi'^X2'^P(xi,X2;t) ,.. , ^ , 1 ^n 

J JdXidX2PixuX2;t) ^W(i,i,tj 



with the product {x1^){x2^). It turns out that the ratio of these quantities depends 
asymptotically on time t, while for any scaling distribution P(xi,X2;t) such a ratio 
would be a constant. In particular, 

(x?)(x^) ■ ^ ' 

Only in the limit n — » this ratio reaches a constant, while for every positive n 
the ratio decays in time. By considering the case n = 1 one sees that the average 
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length, (xi) ~ t^i^-^^"^)/"^ r^ t~-^^^, decays slower than the square root of the average 
area, J{xiX2) ~ t^^^^. This again confirms that the fragment distribution function 
P(xi, X2', t) in the two-dimensional random scission model does not approach a scaling 
form in the long-time limit. However, since all the moments still show the power-law 
behavior we conclude that the model exhibits a multiscaling asymptotic behavior. 

The moments provide an almost complete analytical description of the fragmen- 
tation process. However, the snapshot of the system at the later stages remains 
intriguing (see Fig. 2). This unexpectedly rich pattern arising in such a simple 
process can be viewed as a consequence of the fact that the process is not fully self- 
similar. Instead, the pattern is formed of sets of different scales which are spatially 
interwoven. Fig. 2 also shows that a number of rectangulars have large aspect ratio. 
This qualitative observation is in agreement with the power-law behavior of the n^^ 
moment of aspect ratio, 

{{Xl/X^r) ~ tv^^^-2, (15) 

which is valid for \n\ < 1. Interestingly, this moment diverges when \n\ > 1. The 
aspect ratio appears to be growing in time, in other words, perfect squares have great 
tendency of breaking into long and tin rectangulars. 

If we restrict ourselves to the area distribution function, P{A,t), 

00 00 
P{A,t) = / dxidx25{xiX2 — A)P{xi,X2;t), (16) 



which provides a partial description of our system. We will show that P{A,t) ap- 
proaches the scaling form similar to those found for other one- dimensional fragmen- 
tation systems [y, 0, [1^ . 




Figure 2: Realization of the fragmentation process on a unit square at time t=1000. 



Indeed, the diagonal moments scale in time according to 

s(s + 1 



(17) 



or in other words, the normalized n**^ moments of the area (A")^/" are all proportional 
to t~^. Hence, the area distribution function follows the scaling form 



P{A,t)-t^<^2{At) 
where the scaling function $2(^) satisfies 



6r(s) 



(18) 



(19) 



/o s(s + l) 

Performing the inverse Mellin transforms yields the explicit expression for the scaling 



function fl7 



$2(^) = 6^'dcQ-l)e-^^ 



(20) 



with the hmiting behavior 

^ , , f 62;-2e-^, if z > 1, 

Notice that the scahng solution of Eq. (18) is characterized by the same exponents 
as the scahng solution for the one dimensional random scission model, P{x,t) ~ 
t^$i(xt). However, the scaling functions are different, $1(2) = e~^ (see Eq. (3)) is 
regular everywhere while $2(^) diverges logarithmically near the origin. 

One can consider variations of this model for describing the kinetics of fragmen- 
tation of multidimensional objects. For example, one can change the governing rule 
of the fragmentation events (see Fig. 1) and keep only two rectangulars, say the rect- 
angular in the bottom left corner and in the upper right ones. This rule allows one 
to keep the total length of fragments constant while the total area decays to zero. 
Interestingly, this case has been partially studied in connection with the problem of 
random sequential adsorption of needles [|18| . This model can be treated by apply- 



ing our approach. One should just change in Eq. (4) the factor 4, corresponding to 
creation of four rectangulars, by factor 2. Many results like Eqs. (8), (10), and (12) 
remain the same, with 



Si + S2 , f Si ~ S2^ 



instead of Eq. (9). All the qualitative conclusions also do not change: the model 
exhibits a multiscaling asymptotic behavior and, e. g., 

The area distribution function again scales according to 

P(A,i) ^t^$2(^t). (24) 

Here the scaling function $2 (-2) is given by 

$2(.) = C/'^(1-C)^-V^^ (25) 

Jo C 

where C = r(2v^)/r^(v^) = 2.18482. In the limits of large and small area one has 

r CA-^e-^*, if At > 1, 



^ '^^^ [Dt^lnil/At], ifAt<l, 



(26) 



with D = T{2^/2)/T^{V2) = 2.46432. 

III. Generalization to Higher Dimensions 

We turn now to the general d-dimensional random scission model. The asymptotic 
method presented for the two-dimensional case can be generalized using a simple 
geometric construction. We are interested in obtaining the moment M (s; t) with the 
notation s = (si, . . . , s„). In analogy with the two-dimensional case, we assume a 
power-law behavior M(s; t) ~ t~°^®^ The exponents a satisfy the following difference 
equation 

a(s) + l = a(s+l), (27) 

with the notation 1 = (1, . . . , 1). Meanwhile, the exponents should also reflect the 
hidden conserved integrals, z. e., a(s*) = on the hypersurface s*, HjS* = 2°^. The 
solution to Eq. (^^ with these boundary conditions is given by the formal expression 

a(s) = a(s* + kl) = k. (28) 

This solution clearly satisfies the boundary condition as well as Eq. (28) . Hence the 
problem is reduced to finding roots of the algebraic equation Ilj{sj — k) = 2'^. Since 
this equation is of degree d, a solution is feasible only for d < A. An alternative way 
of viewing the solution is geometric. In Eq. (|2^) s* + kl represents a line along the 
(!,...,!) direction originating at s* and ending at s = s* + kl. The exponent a(s) 
equals the projection of this line on an arbitrary axis, e. g. , on the si axis. Figure 
3 illustrates this construction. 



S2 




Si 



Figure 3: The Geometric Solution. The hypersurface s* satisfies 11, s, = 2 



The main features found for the two-dimensional case such as multiscaling occur 
for higher dimensions as well. As a manifestation of the existence of multiple length 
scales in the system let us consider the ratio of the average volume {V) to the d^^ 
power of the average length, (/). We define the exponent Pd by 



{V)_ 

(0 



-/3d 



(29) 



or equivalently, (3d = l — d (a(2, 1 . . . , 1) + 1). Using the construction of Eq. (^71) , we 
find Pd = 0, 0.1231, 0.1486 for d = 1, 2, 3 respectively, while in the limit rf — »• cxo this 
exponent saturates at 1 — 2 log(3/2) = 0.1891. Note that Pd measures the deviation 
between the asymptotic behavior of the length and the volume. As the dimension 
is increased, this discrepancy becomes more pronounced, and hence, multiscaling is 
stronger in higher dimensions. Another consequence of the same phenomenon is the 
nonuniversal behavior of the various moments of the length distribution. We find 
that the n^^ moment decays asymptotically according to 



(r) n^r'^^°si^+"-m/'^^ 



(30) 



indicating the presence of an infinite number of length scales. 

One can also show that different directions behave independently to a certain 



degree in the limit of infinite dimensions. Specifically, one can show that 

n,x;^) = n,(x?), (31) 



3-^] I '-'-3 \^3 

if Tij = except for a finite number. The average over a finite number of variables 
decouples into a product over single variable averages, while the average over an 
infinite number does decouple. 

For completeness, we present the general dimension results for the diagonal mo- 
ments, M{s, . . . ,s;t), which will be shortly denoted by M{s, t). The governing equa- 
tion for these moments reads 

^^^=[(2/.)'^-l]M(. + l,t). (32) 

We substitute the power-law asymptotic behavior, M{s^ t) ~ A{s)t''°'^^\ into Eq. (32) 
and take into account the boundary conditions a(s = 2) = and A{s = 2) = 1. By 
solving the resulting difference equations we find 

ais) = s-2, A(g) = r-^(g) n r|, _ 2 ■ c^l (33) 

with ( = exp{2'n'i/d). 

In the long-time limit the volume distribution function, P(y,t), approaches the 
scaling form 

P{V,t)-t^<^,iVt), (34) 

with $d(-2) being the inverse Mellon transform of A{s). After a lengthy calculation 
one can find the asymptotic behavior of the volume distribution function: 

' CdV-^e-^\ if Vt > 1, 

Ddt^\n'^-\l/Vt), if\/t<l, 

where Ca = Il'j=l r(2 - 2 ■ C^) and Da = 2'^~i(2'^ - l)/T{d). Thus for all rf > 1 the 
volume distribution function diverges logarithmically in the small- volume limit. 

To summarize, in the d- dimensional random scission model, the volume is charac- 
terized by only one scale, V ~ t^^. However, other geometrical characteristics such 
as the average length and the surface area decay nonuniversaly in time because of 
the existence of an infinite amount of length scales, namely multiscaling. 

One can also consider a varying fragmentation rate and study the case when the 
overall rate depends on the volume as a power-law, i. e., as V^ (the case A = 1 cor- 
responds to the random scission model). When A is positive, this generalization also 
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wt)-^:'^ ,,„,_,,:_ _;_:; (35) 
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Figure 4: The average length of a polygon side, for the random orientation fragmen- 
tation process. Shown are {l{t)) vs. t (diamonds) and a hne of slope —1/2 (solid) for 
reference. 

results in multiscaling of the fragments distribution. The total number of fragments 
N(t) = M{l;t) grows algebraically in time, N{t) ~ t^/'*'. Hence, the case A = is 
a critical one and the number of fragments grows exponentially in time. Finally, for 
A < the shattering transition takes place: the total volume decreases monotonically 
and the total number of fragments reaches infinity within an infinitesimally small time 
interval. Moreover, a finite fraction of the volume breaks into zero-volume rectangu- 
lars. This phenomenon is well known in the context of one- dimensional fragmentation 
[p!0| , pT| and has been examined in the context of two-dimensional fragmentation with 



length conservation in a very recent study [jT9 . 

IV. Isotropic Fragmentation 

Intrigued by the rich kinetics of the rectangular fragmentation problem, we also 
investigated numerically an isotropic fragmentation process. In situations such as 
shattering of a thin glass plate or in membrane crumpling, the fragments are polygons 
with a varying number of sides. Hence, we introduce a process where a randomly 
oriented crack appears with a uniform rate at a random point of the surface and 
propagates with infinite speed until it meets existing cracks. The original model 
can be viewed as deposition of such perfectly oriented "cross" shaped cracks. The 
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overall fragmentation rate in both processes equals the volume of the fragment. For 
randomly oriented fragmentation, each fragmentation event creates an additional 
polygon and hence the total number of polygons grows linearly in time according to 
N{t) = M(l, 1; t) = 1 + 1. The average volume thus scales as 1/N{t) or A ~ t'^. 

A Monte-Carlo simulation study of isotropic fragmentation process on a unit 
square suggests that unlike for oriented fragmentation, only one length scale exists 
in the isotropic problem. The average length of a polygon side is plotted in Fig. 4 
and appears to decay as t^^/^. Therefore, the length follows the same asymptotic 
behavior as does the square root of of the average area. A snapshot of a realization of 
the system at time t = 1000 is shown in Fig. 5. This interesting picture suggests that 
it might prove insightful to investigate various structure properties of the system such 
as the area distribution function and the side number probabilities of the polygons. 

In conclusion, we have studied two fragmentation processes in spatial dimensions 
larger than one. For oriented fragmentation, where the fragments are always rectan- 
gular, multiscaling is found in the long-time limit. Specifically, the length distribution 
function has moments that scale algebraically in time with an infinite number of in- 
dependent length scales, while the area distribution function is characterized by a 
single length scale. The area distribution function also exhibits a weak logarithmic 
singularity near the origin. Multiscaling appears to depend strongly on the geometric 
nature of the process. For isotropic fragmentation, a single length scale describes the 
decay of the length as well as the area. 
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Figure Captions 

1. Illustration of the fragmentation process in two-dimensions. 

2. Realization of the fragmentation process on a unit square at time t = 1000. 

3. The geometric solution, the hypersurface s* satisfies HjSj = 2'^. 

4. The average length of a polygon side, for the random orientation fragmen- 
tation process. Shown are (/(t)) vs. t (diamonds) and a line of slope —1/2 
for reference. 

5. Realization of the random orientation fragmentation on a unit square at 
time t = 1000. 
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